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Variations in Atmospheric Angular Momentum 
and the Length of Day 

Richard D. Rosen and David A. Salstein 

Twice-daily global analyses of the zonal winds prepared by the 
U.S. National Meteorological Center for the period from 1 January 1976 
through 31 December 1981 have been used to construct a time series of 
the atmosphere's angular momentum (M) . Three-day means of this quantity 
for the first five years of the data set are spectrally analyzed and 
compared with independent observations of changes in the solid earth's 
rate of rotation. The most prominent feature in the M time series is 
an annual signal, and this is mirrored by similar behavior in changes 
in the length of the day (l.o.d.). Other noteworthy fluctuations in M 
and l.o.d. on time scales of less than a year also occur, however, and 
these too agree well with each other. In this regard, a near 50 day 
fluctuation in both time series, also noted in other studies, is espec- 
cially visible. We conclude on the basis of these comparisons that 
indeed most of the variability in our time series for M is real and 
that, as others have suggested, the atmosphere plays a dominant role 
in forcing changes in l.o.d. on time scales of about a year and less. 

We divide the mean meridional plane into a number of regions to 
investigate the contributions made by different parts of the atmosphere 
to the changes found in the global statistic M. Our approach is to 
calculate empirical orthogonal functions that define the major spatial 
modes of variation among the regional time series of angular momentum. 
As expected, we find the annual cycle in M is associated with seasonal 



changes in the major jet streams of the two hemispheres, with the larger 
of these changes occurring in the Northern Hemisphere. We also find 
that all the statistically significant spatial modes of non-seasonal 
variation depict behavior organized over large horizontal scales with 
a high degree of vertical stratification. Non-seasonal changes in M 
emerge from the imbalances among these regional variations in momentum, 
with the largest contributions on time scales less than a year coming 
from the subtropical portions of each hemisphere. 
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1 . 


Introduction 


A fundamental measure of the dynamic state of the atmosphere is 
its angular momentum about the polar axis relative to the earth. Indeed, 
study of the maintenance of the atraosph re’s angular momentum against 
dissipation has long been at the forefront of n*eteorological research 
(see Lorenz, 1967, for a review of this history). In an early, ar.d now 
classic, paper on the topic, Starr (1948) also noted that M there is no 
reason to expect that the partition of angular momentum [between the 
earth and the atmosphere] should remain constant when seasonal and other 
short time-intervals are considered”. Since then, numerous studies have 
been undertaken to quantify more precisely the relationship between 
changes in atmospheric momentum and the rotation rate of the solid 
eartn on a variety of time scales. 

Actually, the rotation of the solid earth is affected by a wide 
range of geophysical phenomena besides those related to the atmosphere. 
These include the effects of tides, earthquakes, core-mantle coupling 
(Yoder et al., 1981; Lambeck and Cazenave, 1977; Hide, 1977) and per- 
haps the melting of polar ice (Etkins and Epstein, 1982). On periods of 
less than about a year, however, the atmospheric excitation of changes 
in the earth’s rotation rate appears to dominate (except at the fort- 
nightly and monthly solid body tidal periods). Indeed, the atmosphere 
seems to play a considerable role in this regard out to periods of aoout 
four years (Lambeck and Hopgood, 1981, 1982), wh »xh has made more pre- 
cise determinations of the atmospheric contribution of great interest 
to geodesists and geophysicists. 

Most of our previous knowledge concerning temporal changes in the 
atmosphere’s angular momentum has been derived from monthly mean zonal 
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wind data at the network of upper-air sounding stations. These stations, 
though, are irregularly spaced and for the most part are located over 
the land north of 20° S, so that their representativeness has been a 
concern. Within the last decade, however, several large meteorological 
centers have introduced global analyses of the wind field that routinely 
assimilate data not only from the rawinsonde stations, but also from 
commercial aircraft and satellites which provide important new informa- 
tion in the station-sparse regions. In addition, these global analyses 
are convenient to use, since they provide wind values directly on a 
latitude-longitude grid. In particular, the computation of atmospheric 
angular momentum on a daily basis is straightforward. Recently, Hide 
et al. (1980) demonstrated the potential that now exists for studying 
high frequency changes in atmospheric angular momentum by using the 
global grid point analyses produced by the British Meteorological Office 
and by the b.S. National Meteorological Center (NMC) for the same four- 
month period. We extend the work of Hide et al. here by using six years 
(1976-81) of NMC twice-daily global analyses to create and study a leng- 
thy time series of high temporal resolution angular momentum values. We 
also compare changes in these atmospheric values to independently deter- 
mined changes in the rotation rate of the solid earth. Finally, we ex- 
amine the atmospheric data in more detail to determine the time and 
space scales on which variations in momentum occur within the atmosphere 
and which regions are contributing most to the changes we find In the 
global integral. We begin in the next section, though, with a descrip- 
tion of the data and techniques used to derive our t*me series of 


momentum values. 



ORIGINAL PASS IS 
m POOR QUALITY 


2 . Derivation of Atmospheric Angular Momentum Values 

NMC introduced its first operational glcoal analysis scheme in 
September 1974. Since then, numerous modifications to the scheme have 
been introduced (a detailed description of a recent version is provided 
by McPherson et al. , 1979), but from the start the output has included 
twice-daily (0000 and 1200 GMT) values of the zonal wind analyzed over 
a grid with points spaced every 2.5° in both latitude and longitude, at 
each of 12 pressure levels in the vertical (1000,850,700,500,400,300, 
250,200,150,100,70 and 50 mb). NMC has archived certain general circu- 
lation statistics derived from these grid point analyses beginning 1 Jan- 
uary 1976 (Miller et al. , 1975). Included in this archive are zonally 
averaged values of the eastward component of the wind at each latitude 
and level, [u]. To derive estimates of the angular momentum (M) of the 
atmosphere about the polar axis, relative to an earth-fixed frame, these 
[u] values have been integrated, with the appropriate weighting, over 
latitude ($) and pressure (p) : 


M 


27ra 3 

g 


100 mb —it / 2 

[u] cos 2 <j> d$ dp 

1000 mb it / 2 


( 1 ) 


where a is the mean radius of the earth (6.37 x 10^ m) and £ is the 

-2 

acceleration due to gravity (9.81 ms ). Equation (1) is based on the 
assumption that the atmosphere is in hydrostatic equilibrium and ignores 
the variation with altitude and latitude of the distance of a parcel of 
air from the center of the earth. We have not used any of the wind data 
at the two levels above 100 mb where the reliability of the NMC analysis 
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is questionable. The trapezoidal rule was used in evaluating (1) 
numerically. 

Occasional brief gaps exist in the NMC archive and no attempt has 
been made to fill these in. However, there are also two relatively 
lengthy gaps, one from 1 April through 15 August 1977 and one from 
1 February through 8 April 1981. For these periods, we obtained copies 
of the original NMC grid point analyses (for 0000 GMT only) from the 
National Center for Atmospheric Research, and we used these to fill in 
the missing values for [u] and M. 

Errors in our estimates of global momentum values arise through 
the approximations involved in deriving (1) as well as through inaccur- 
acies in the NMC wind analyses. With regard to the former source of 
error, the most serious aspect is undoubtedly the neglect of the upper 
atmosphere at pressures less than 100 mb. We have examined some calcu- 
lations made with and without data from the 70 and 50 mb levels, and 
based on these we conclude that neglecting the upper 10% of the atmo- 
sphere incurs a systematic underestimate in the mean level of M of about 
10% or less, but it has a much smaller impact on our determination of 
most short-term changes in M. (On the other hand, some variations in M 
such as a semi-annual or a quasi-biennial cycle may well be affected by 
the lack of data from much of the stratosphere and above.) 

An error also results from our treatment of the lower boundary in 
(1). In the NMC analysis, the lower boundary is taken to be smooth and 
fixed at a constant pressure of 1000 mb, regardless of where this level 
lies relative to the earth’s surface at a particular grid point on a 
given day. When the surface of the earth is actually at pressures less 
than 1000 mb, the NMC analysis will nevertheless place non-zero winds 
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beneath the topography. These spurious winds have been Included in our 
[u] data, but a test in which we removed them from a short period of 
data revealed that their presence alters M by no more than about IX. 

On the other hand, when the 1000 mb level happens to lie above the 
earth's surface, our approach neglects the contribution to M made by 
the part of the atmosphere at pressures greater than 1000 mb. Because, 
however, the average surface pressure over the globe is around 985 mb 
(Trenberth, 1981a), thi^ error is small. 

As we noted above, another source of error in our values of M 
results from inaccuracies in the basic NMC grid point wind analyses 
themselves. Although NMC incorporates measurements from an extensive 
observing network, large data gaps still exist, particularly over the 
oceans. Therefore, the analyses over such areas are bound to be im- 
perfect. Nonetheless, our examination of an early version of the NMC 
global analysis over the Northern Hemisphere did not reveal any obvious 
defects in the zonal wind field (Rosen and Salstein, 1980). Some aspects 
of the analysis error in a more recent version of the NMC analysis are 
presented by McPherson et al. (1979), but it is difficult to estimate 
from their discussion what the magnitude of the error in M is likely to 
be. Oort (1978) examined the adequacy of the rawinsonde network for 
determining circulation statistics and found that the most significant 
error in M was on the order of 5% and caused by the presence of spatial 
gaps in this network. Although Oort's study was not based on the NMC 
analysis and did not account for the effect of satellite observations, 
his result does provide some basis for estimating the size of the random 
errors in our M values. The comparisons between the NMC analysis and 
that of the British Meteorological Office made by Hide et al. (1980) 




are also revealing in this respect. Differences as large as 10% 
between the two analyses did occasionally occur, but much of this 
appeared to be systematic in nature. 

Occasional procedural changes at NMC during the six-year period 
introduce an unknown degree of heterogeneity to our data. Major 
changes occurred ir September 1978 and May 1980 when new analysis 
methods were introduced (again see McPherson et al., 1 97 9 » also Sela, 
1980), but other important modifications have also been made that were 
not always well documented. In addition, ther^ have been continuing 
changes in the satellite observing systems used by the NMC analysis. 

For the most part, though, we feel that these sorts of changes have 
little impact on day-to-day variations in our M values. 

In summary, we believe that we may be systematically underestim- 
ating the mean value of M by about 10%. Random errors that affect our 
estimates of most short-term changes in M are probably less than 10%. 

3. Results for Atmospheric Angular Momentum 

All available values of M for each of the calendar 'ars 1976-1981 
are plotted in Fi^s. l(a)-(f). In addition to the global atmospheric 
momentum values, time series are also plotted for the angular momentum 
of the atmosphere above the Northern Hemisphere (NH) anc the Southern 
Hemisphere (SH) . (The scale on the right in the figure gives the equiv- 
alent measure of M in terms of inferred changes in the length of day; 
more on this in section A.) 

Perhaps the most striking feature displayed in Fig. 1 is the 
annual signal in the global M values. Largest values occur in December 
or January, and the smallest in July or August. This pattern occurs 
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because the momentum time series of the two hemispheres are far from 
being mirror images. Although the annual wave in each hemisphere peaks 
in its respective winter season, its amplitude is clearly much larger in 
the NH. In the SH, the peak in the time series is generally quite broad 
and values in the SH winder are typically only 25% larger than in the 
summer. Presumably, this contrast between the two hemispheres is related 
to their very different distributions of land and ocean. Our results for 
the SH appear to be supported by those derived by Swanson and Trenberth 
(1981a) from the SH geopotential height analyses of the Australian Wea- 
ther Bureau. They find that the SH winter westerly geostrophic jet is 
only 25% stronger than the summer jet at upper levels, whereas at 500 mb 
and below the maximum winds actually occur during the SH summer (albeit 
at more northerly latitudes than in the winter, thereby diminishing their 
effect on M). A final note regarding the hemispheric differences in Fig. 
1 is that instances occur during which the relative momentum of the atmo- 
sphere above the NH is actually negative; this never happens over the SH. 

The dominance of the annual signal in the time series for global 
M is clearly evident when we compress our data into a single plot. This 
is done in Fig. 2, wherein present all our M values for the five-year 
period 1 January 1976 through 31 December 1980. ^ One feature of the 
annua] signal made especially visible in Fig. 2 is its departure from a 
pure sinusoidal form, with its peaks tending to be broader than its val- 
leys. It is clear, too, that interannual variations in M exist. An 
obvious example is afforded by the anomalously low M values during the 
first part of 1976. Inspection of Fig. 1 reveals this is due mainly to 

^Because we acquired the data for 1981 just recently, they have not been 
included in Fig. 2 nor in any of the subsequent analyses presented here. 
The data for 1976-80 are tabulated in Rosen and Salstein (1981). 
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behavior in the SH, although NH momentum values are also somewhat de- 
pressed at this time. Since observations are particularly sparse over 
much of ‘.he SH, one may well question the reality of this anomaly. On 
the other hand, Krueger (1982) has analyzed monthly mean values of 
global M back to early 1975 and finds that the behavior of M during the 

period from late 1975 through 1976 resembles that of some other atmo- 

2 

spheric indices including a measure of the Southern Oscillation. Also, 
using analyses by the Australian Weather Bureau, Swanson and Trenberth 
(1981b) find that the geostrophic westerly zonal wind component was 
anomalously weak during early 1976 throughout the depth of the SH sud- 
tropics. Although there is some indication (Anderson, 1982) that the 
NMC wind analysis in higher SH latitudes may have been faulty dt g 
this period, we have chosen in the face of inconclusive evidence simply 
to accept the M data for 1976 as they stand here. 

Higher frequency variations are, of course, also present in the M 
time series. Returning to Fig. 1, for example, a small diurnal component 
is visible when most of the curve is contrasted with a period when only 
once-daily data were available such as April through mid-August 1977. 

A fairly large fluctuation of about 50 days may also be seen at times, 
for example during 1981 when it is particularly strong. We wiJ 1 have 
more to say about this 50-day fluctuation in later sections, but first 
let us consider more generally a time series of M from which the large 
annual signal has been removed. To construct such a series, we began by 
forming the average of all available 0000 GMT values of M in each succes- 
sive 3-day period for 1976-80. Aside from filtering out contributions 

^Utilizing monthly mean winds for May 1963-April 1973, Stefanick (1982) 
also finds a relationship between M and the Southern Oscillation. 
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by very high frequencies (according to Table 1, fluctuations with per- 
iods of less than 3 days contain only about 2% of the total variance in 
the full M time series), this procedure also provides a convenient means 
for bridging gaps in the o^ginal time series. Thus, although 10% of 
the once-daily data were missing during the 5-year (1827 day) period, 
only five of the ( 09 3-day periods had no data (M values for these were 
assigned by interpolating surrounding 3-day means). 

To define the annual signal, we averaged the corresponding 3-day 
means from each of the five years. This signal was then subtracted 
from the original 3-day mean values to yield the result in Fig. 3. The 
variance remaining in this residual data set is about 15% of that con- 
tained in the original time series of M displayed in Fig. 2. Much of 
the variability in Fig. 3 appears to be a.ssociated with the near 50-day 
fluctuation noted earlier. Clearly, its phase and amplitude change con- 
siderably throughout the record. 

A quantitative assessment of the frequencies present in Fig. 3 
may be obtained, of course, through use of soectral analysis techniques. 
Because we wished to place confidence limits on our estimate of the spec- 
tra and because our data are limited to a five year period, we decided to 
adopt an analysis technique that would focus on fluctuations with periods 
shorter than about 100 days, roughly 1/18 of the length of our record. 
Our approach follows that outlined by Welch (1967). first organized 

our time series of 3-day mean values of M for 1976-80 into eight over- 
lapping blocks spanning most of these five years. Each block consists 
of 128 consecutive values of the 3-day means multiplied by a weighting 
function suggested by Welch which sets the endpoints of each block to 
zero. A fast Fo ur ier transform method was then applied to determine the 
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Mean, standard deviation and variance of global M during 1976-80 (00 GMT only) 
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power spectrum in each block. Lastly, the resulting eight spectra were 
averaged to produce a final result. 

By organizing our data into smaller blocks, the possibility is 
lessened that * aifts in the phase of a particular frequence auring the 
course of the entire record will act to reduce the power perceived at 
that frequency. Another advantage of the approach followed ‘.ere is that 
it permits the statistical significance of peaks in the spectrum to be 
readily assessed. We do so y first making a power-law best fit to the 
spectrum (except for its highest frequencies), in concert with the known 
(Ward and Shapiro, 1961) red-noise nature of most meteorological time 
series. We then obtain confidence limits about this best-fit curve by 
applying Welch 1 s model of the spectrum expected from random noise. 

We present the final results of our spectral decomposition of M 
for oscillations with periods of from 6 to 128 days in Fig. 4. Spectra 
of atmospheric momentum fluctuations on periods of two months and longer 
were presented by Lambeck and Hopgood (1981) based on monthly mean values 
of M. Our data, of course, can reveal higher frequency components, and 
as the figure demonstrates several peaks in the 6-128 day range, such as 
the prominent one near 30 days, do approach or exceed our 93% level of 
statistical significance. Of particular note is the rather broad peak 
centered at about 50 days, confirming our earlier qualitative remarks 
about the presence of this signal. (The strength of the 50-day peak 
would likely have been even larger had we included the 1981 data in our 

analysis, given the obviously strong presence of this signal in Fig. If.) 
We should remark, of course, that passing a test of statistical signif- 
icance does not guarantee that a particular fluctuation results from an 
underlying physical mechanism. Additional, support for the reality of 
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the 50-day fluctuation, however, is provided by studying independent 
measures of the rotation rate of the solid earth. More about this 
matter is discussed in the next section. 


4. Relationship Between Changes in the Momenta of the Atmosphere 

and Solid Earth 

As we discussed at the outset, there have been a number of studies 
relating changes in the atmosphere’s momentum to changes in the rotation 
rate of the solid earth. The simplest model relating these two is to 
treat the (rigid) earth and atmosphere as a closed dynamical system and 
assume that a change AM in the angular momentum of the atmosphere is 
accompanied by an equal, but opposite, change in the angular momentum 
of the earth. The consequent change in the angular velocity u) of the 
surface of the earth is given by: 



( 2 ) 


where I is the principal (axial) moment of inertia of that portion of 
the earth that responds to the change in M on the time scale being con- 
sidered. Here we are mostly interested in time scales on the order of 
one year or less, and fluctuations on these periods are thought to affect 
only the earth’s crust and mantle (i.e., its shell) and not its core 

(Hide et al. , 1980). Following Langley et al. (1981b), therefore, we 

37 2 

have set I * I , n * 7.04 x 10“ kg m . 

shell ° 

It has become customary to speak of changes in the solid earth’s 
rotation rate not in terms of A.o but rather in terms of changes in the 
length of day (Al.o.d.). The two are simply related by: 
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Al.o.d. m -Au) ^ -Am 
1 • o . d * m + Aw w 

Combining (2) and (3) and setting l.o.d. ** 86400 s, w ■ 7.29x10 ^ s ^ 
and I as above, we find the following linear relation between AM and 
the Al.o.d. associated with it: 


Al.o.d. - - - AM = 1.68x 10* 29 AM (4) 

wl 

2 -1 

where Al.o.d. is in units of seconds and AM is in kg m s . Noting 

r* f 

from Table 1 that the mean value of M during 1976-80 was VI. 4 x 10 kg 
2 -1 

m s , we may infer from (4) that were this amount to be transferred 
entirely to the earth’s shell the l.o.d. would change by 2.3 ms. 

For convenience, we have chosen to reckon AM as changes in M from 
a hypothetical base state in which M = 0. Using this definition in (4) 
allows us to equate the time series of M values directly in terms of 
their implied changes in l.o.d. from this base state, and we have done 
so in the scale drawn along the right-hand ordinate of Figs. 1-3. The 
question naturally arises then as to how closely these inferred values 
of Al.o.d. compare with those determined directly through astronomical 
measurements. 

The classical technique for measuring changes in the earth's 
rotation is based on observing the transit times of stars at a network 
of astronomical observatories. This approach and the results derived 
from it have been reviewed by Larabeck (1980). Recently, two other 
space-based techniques have begun to provide values of Al.o.d. One of 
these is very long baseline interferometry (VLBI), in which different 
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antennas receive radio emissions from the same distant quasar (Robert- 
son and Carter, 1982). The other technique involves measuring the time 
required by a laser beam to travel from the earth to a satellite and 
reflect back. The satellite can either be a specially designed arti- 
ficial one, as in the case of the Lageos satellite (Smith et al., 1979; 

Schutz et al., 1981), or the Moon, which was equipped with mirrors 
placed by the astronauts (Mulholland, 1980; Langley et al. , 1981a). 

Values of Al.o.d. at 5-day intervals for the 1976-80 period determined 
by lunar laser ranging (LLR) have been kindly provided to us by R.B. 

Langley and colleagues at M.I.T. and, therefore, serve as a convenient 

data set for us to use. The fortnightly and monthly tidal terms have 

been removed from these data, and the resulting values are plotted in 

Fig. 5 along with cur own time series of Al.o.d. inferred from M through (4). 

We should first point out that the zero level on the scale for 
each curve in Fig. 5 is arbitrary. The curve using the M data is, as 
we noted above, reckoned from the hypothetical base state in which M ■ 0; 
the curve using the LLR data measures the length of day in excess of that 
of the mean solar day during the nineteenth century, which by convention 
is 86400 s (Hide et al., 1980). The mean levels of the two curves in the 
figure do happen to be close, but of course we should not conclude from 
this that M was zero on the average during the nineteenth century. 

Rather, it is clear that torques other than those imparted by the atmo- 
sphere must be important in affecting the rotation rate of the shell on 
these decade and century-long time scales. Even for the five year period 
shown in Fig. 5, the LLR curve shows evidence of a trend that is not 
apparent in the M curve, again implying that (4) oversimplifies the 
situation that exists on time scales beyond about a year. 
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On time scales shorter than a year or so, however, the two curves 
in Fig. 5 agree quite well. Certainly they both contain similar annual 
signals. The resemblance between the two curves at higher frequencies 
may be more strikingly displayed if the mean levels and annual and semi- 
annual variations are removed from each and the residuals smoothed, as 

3 

Langley et al. (1981b) do for each of the years 1976-79 . We reproduce 
their result here as our Fig. 6. The most prominent feature in this 
figure is the close correspondence in the amplitude and phase of the 
two curves on periods of about 50 days. The agreement is particularly 
good during 1979, the year of the Global Weather Experiment in which 
more extensive observations may have resulted in improved NMC wind 

analyses. Interestingly, the other astronomical measures of Al.o.d. 
mentioned before also exhibited the nearly same ^50-day fluctuation 

during 1979 (Feissel and Gambis, 1980). We gain confidence, therefore, 
in the reality of this fluctuation in our M time series. In fact, we 
can more generally conclude on the basis of the comparisons here (and 
in Feissel et al. , 1982, between our M data and the classical technique 
for measuring Al.o.d.) that most of the variability in our M data is 
realistic, at least on time scales from about a week to a year where 
support is offered by the independent data sets. We may also conclude 
from Fig. 6, as do Langley et al. (1981b), that the differences between 
the two curves lie within their combined uncertainties and that contri- 
butions from non-atmospheric torques on the time scales displayed can 

4 

probably be discounted. 

guassian- shaped low-pass filter with a window width of about 8 days 
was used in the smoothing. 

^Recall, though, that the large effects associated with tidal induced 
changes in the earth's moment of inertia on fortnightly and monthly 
periods have been removed from the observed (LLR) changes in l.o.d. 
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Finally, we should note that the approach followed here does not 
shed any new light on the manner in which momentum is exchanged between 
the atmosphere and solid earth. We have simply demonstrated that such 
an exchange must be taking place. Swinbank (1981), however, has recently 
examined the roles played by the two major torques (surface friction and 
mountain) that must be responsible for the momentum exchange. Using 
daily data for the two special observing periods of the Global Weather 
Experiment (a total of about four months during 1979) and a suitable 
model for the frictional torque, he shows that most of the variability 
in the rate of change in M is due to the mount in torque. Friction, 
although of comparable magnitude to the mountain torque, is however 
much steadier, and it does not correlate at all well with the time 
derivative of M. Naturally, it would be desirable to extend Swinbank' s 
study to other periods and to other models of these torques. 

5. Variations of Angular Momentum within the Atmosphere 

a. Variations in latitude-height cross sectional boxes 
Thus far we have dealt with the atmosphere's momentum on either 
global or hemispheric spatial scales. To understand more about the 
reasons for the fluctuations observed in these large-scale statistics, 
however, it is necessary to examine the atmosphere on smaller scales. 

For example, the origin of the 'vSO-day variation in M cannot be explained 
on the basis of these global values alone. On the other hand. Madden and 
Julian (1971, 1972) have reported that such a periodicity existed in sev- 
eral tropical station upper-air wind reports during 1957-67,^ thereby 

^It was on the basis of this work by Madden and Julian, in fact, that 
Lambeck and Carenave (1974) first suggested there might be atmospheric 
forcing of a v*50-day fluctuation in l.o.d. 
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providing a regional focus for furthe’ *tudy of this particular phenom- 
enon, We are, of course, interested in studying initially the entire 
broad range of fluctuations present in our atmospheric data, and so we 
have not specially isolated any one periodicity for detailed ccrutiny 
here. 

Our first step in studying the structure and behavior of the field 
of momentum within the atmosphere was to compress our zonal mean wind 
data set. With values of [u] available a r 73 latitudes and 10 pressure 
levels, we could conceivably deal with 730 different time series. This 
amount of data, however, would have proven too cumbersome for the anal- 
ysis approach we planned to US e, an ^ so we took advantage of the known 

spatial cohesiveness of large-scale atmospheric motions to congregate 
the 730 grid points into 27 equal mass "boxes." These boxes were formed 
by dividing the atmosphere into three equal mass layers in the vertical 
(1000-700 mb, 700-400 mb, 400-100 mb representing the low, middle and 
upper troposphere, respectively) and nine equal area belts in the hori- 
zontal with boundaries at 6.4°, 19,5°, 33.7°, 51.1° and 90° in each hem- 
isphere. For simplicity, we will refer to the zones delimited by these 
latitude circles in the following way: 6.4°S - 6.4°N, equatorial; 

6.4° - 19.5°, tropical; 19-5° -33.7°, subtropical; 33.7° - 51.1°, mid- 
latitude; and 51.1° - 90°, high latitude. We recognize the possibility 
that this scheme may stratify our data too severely to investigate prop- 
erly the structure of the momentum field. Preliminary calculations with 
the full 730 grid point data set by us and work b> Anderson (1982) with 
these same data suggest, however, this is generally not the case. 
Finally, we note that the zones and layers we are using do correspond 
roughly to the broad climatological regions of the general circulation. 
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We next computed the angular momentum of the atmosphere within each 


box, B, viz: 


| cos2 (5) 

B 

where dA is an elemental area of the box in the pressure- latitude plane. 

We will refer to ra^ as the momentum density. To evaluate (5) numerically, 

we applied the trapezoidal rule in the pressure and latitude, coordinates 

piecewise between the 730 data points, being careful to ensure that * M 

27 boxes 

each day. 

We have chosen here to work entirely with 3-day means of the once- 
daily (00 GMT) values, which we formed in the same manner described 
in section 3 to construct 3-day means of M. A plot of the 5-year mean 
momentum density field based on these data for 1976-80 is given in Fig. 

7. The figure reflects the well-known main features of the general cir- 
culation: low-level easterlies in the equatorial and tropical regions, 

and zonal mean westerly jets farther poleward in the upper troposphere 
of each hemisphere. The variance of the m fi time series in each box for 
1976-80 is displayed in Fig. 8. The maxima in this field are clearly 
associated with the zonal mean jets. Moreover, it is evident that al- 
though m^ is larger in connection with the SH jet, it is more variable 
in the region of the NH jet. 

We wish to examine the spatial and temporal relationships that link 
the behavior of the m time series in the different boxes and give rise 

D 

to the variance values in Fig. 8. To this end, we have adopted an anal- 
ysis approach, known as empirical orthogonal function (EOF) or principal 
component analysis, that yields modes which can identify patterns of 
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large-scale variability, A major advantage this approach has over the 
more conventional one of expanding data in terms of a pre-chosen set of 
analytic functions is that the EOFs explain the variance in the data in 
the most efficient manner possible. We will not detail the basis of the 
approach here (Lorenz, 1956 and Kutzbach, 1967 are excellent references 
for this) , but we will briefly outline the manner in which we have ap- 
plied it. First we formed a 27 x 27 covariance matrix by removing the 
mean value from each m^ time series and then multiplying the resulting 
nig anomaly series in each box by itself and by that in every other box. 
This matrix is then diagonalized and its eigenvalues and eigenvectors 
(the EOFs) found. The anomaly data are then projected onto each eigen- 
vector to determine the time series associated with it. The amount of 
the variance in the original data explained by each eigenvector is equal 
to the variance in the vector* s associated time series, but it is more 
usual to compute the percent of the total variance explained by each 
vector by simply taking the ratio of its eigenvalue to the sum of all 
the eigenvalues. 

Fig. 9 presents the first mode of the variability in the momentum 
density field for 1976-80 and its associated time series. It is clear 
from the time series that this mode is mostly depicting the annual sig- 
nal in the field. The vector, which explains 81.5% of the total vari- 
ance in the data, displays two main centers of action, one in the sub- 
tropical upper troposphere of each hemisphere.^ This result and the 
fact that the NH center is the larger may, of course, have been antici- 
pated on the basis of the picture in Fig. 8 of the m^ variance. Note 

6 Each EOF mode has been normalized so that the sum of the squares of 
the 27 non-dimensional box values is unity. 
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finally with regard to Fig, 9 that the two centers have opposite signs , 
so that when the momentum is larger than average in the NH (as in the 
NH winter when the positive value in the NH center is multiplied by a 
positive value in the time series), it is smaller than average in the 
SH (since the positive value in the time series now multiplies a nega- 
tive value in the SH center) and vice versa. 

The second EOF is presented in Fig. 10 along with its time series 
and a spectrum of the time series determined in the same manner used to 
derive the global M spectrum. This mode explains 5.7% of the variance 
in the m^ data, or nearly 1/3 of the variability left unexplained by the 
first mode. Unlike the first mode, the major contribution to mode 2 is 
made by the SH tropical and subtropical upper-troposphere regions, 
with thest- acting in phase with their counterparts in the NH. Many of 
the statistically significant periodicities evident in the global M 
spectrum (Fig. 4), including the one near 50 days, reappear here in con- 
nection with this EOF. 



further, we have subjected the ro residual data to the same sort of EOF 

a 

analysis we applied above to the total (i.e., including annual signal) 

m u data, 
o 

Although 26 independent EOFs can be derived in our analysis, many 
of these can be expected to depict mostly noise in the data that happens 
by chance to be correlated among some of the 27 boxes during the (finite) 
stud/ period. We. of course, are interested in studying only those EOFa 
that are likely to owe their existence to the presence of some under- 
lying physical mechanism. Recently, we (Salstein et al., 1982) and 
Overland and Preisendorfer (1982) have stressed the importance of deter- 
mining the statistical significance of individual EOFs to separate those 
that might be meaningful from those that are not. Our approach here is 
to construct for each box a time series of normally distributed random 
numbers whose mean and variance are equal to those of the real resid- 
ual data. These 27 independently derived series of random nisnbers are 
then submitted to an EOF analysis yielding an eigenvalue corresponding 
to each mode. After repeating this experiment a number of times with 
different sets of random numbers, we noticed that the collection of 
eigenvalues corresponding to each mode number tended to cluster about 
the same value. We decided, therefore, simply to select the largest 
eigenvalue obtained for a mode from 10 experiments as the level of sta- 
tistical significance required for that mode number. As Fig. 12 illus- 
trates, the first six modes of the m D data with the annual signal removed 

D 

passed this test for significance. Beyond mode 6, however, the eigen- 
vectors explain no more of the variance in the data than would vectors 
obtained from a random set of numbers. We, therefore, discuss only modes 
1 through 6 here. 
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These pictures of the modes of non-seasonal variability in the 
momentum density field are contained in Fig. 13. The time series 
associated with each of the modes is presented in Fig. 14, and a spec- 
trum of each time series in Fig. 15. Not surprisingly, mode 1 in Fig. 

13 resembles the picture in Fig. 10 of the second EOF for before the 

annual signal was removed. The resemblance is not perfect, however 
(r specially when their respective time series are compared). The spec- 
tra are quite similar, though, indicating the expected result that the 
annual signal as we define it contains little power at higher frequen- 
cies, We prefer to view the vector in Fig. 13 as the better for depict- 
ing 'the major mode of non-seasonal variations in m^. The maxima in this 
mode are located in the same regions as the maxima in the variance field 
of the filtered data (Fig. 11), although in the EOF the magnitudes of 
the SH centers are larger than the one in the NH. The large values of 
mode 1 (and the variance) in the SH mid-latitudes appear to be related, 
at least in part, to the anomalous behavior in M during 1976. In light 
of our earlier discussion about the uncertainty in the reality of this 
anomaly, we suggest caution in relying too strongly on our result,*; for 
the SH mid-latitude boxes. Finally, with regard to the spectrum for 
mode 1, we should note especially the presence of the significant peak 
near 50 days. Unfortunately, the analysis approach we have chosen here 
is inadequate for isolating further which subset of the boxes is most 
responsible for this particular periodicity* The work by Madden and 
Julian referenced earlier and the recent study by Anderson and Rosen 
(1982), however, strongly suggest the large values in the tropical boxes 
in Fig. 13a are related to this phenomenon. 


The sotrewhat muted role in mode 1 played by the Nil sub-tropics, 
despite the maximum there in the m s variance, appears to be largely 

O 

compensated for by this region's role in mode 2 (Fig. 13b). The value 
reached by the l subtropical, upper- troposphere box in this mode is, 
in fact, the largest to be found in any of the fields in Fig. 13. The 
striking feature about the spectrum of the mode 2 time series is its 
large peak around 30 days. Given the major role of the NH subtropics 
in this mode (the three boxes in this belt account for 64% of the var- 
iance in the mode) , we can infer this region is responsible in large 
measure for this periodicity. 7 Interestingly, the next most important 
region in the mode 2 vector is the SH mid-latitudes P which whe same 
sign as the NH subtropics. Trenberth (1981b) has demonstrated that this 
area does contain a considerable amount of variance on time scales of 
8-64 days in such fields as 500 mb geopotential height and kinetic en- 
ergy. Our results for this region here for mode 2 and earlier for mode 
1 are not inconsistent with this. 

Although modes 3 through 6 explain successively smaller amounts of 
the variance in m^, they nonetheless depict large-scale cohesive patterns 
of behavior. In mode 3, the NH mid-latitudes are highlighted for the 
first time. The very significant peak in the mode 3 time series spec- 
trum near 13-14 days is due mostly to behavior in this region. This 
periodicity corresponds closely to that found by Miller (1974) in time 
series of various measures of the troposphere's energy cycle north of 
20° N, and attributed by Miller to lunar tidal forcing mechanisms and/or 
amplitude vacillation of tropospheric wave patterns. 

7 An independent check of the mg time series spectrum for the NH sub- 
tropical, upper-level box confirms the presence of a significant peak 
near 30 days. 
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The mode 4 vector is noteworthy for its remarkable degree of sym- 
metry about the equator and striking banded structure. In mode 5, the 
maximum in the SH high latitudes is out of phase with the rest of the 
hemisphere. uterestingly , the spectrum of this mode is distinguished 
by its lack of strong peaks. In mode 6, the SH high latitudes are again 
at a maximum, but this time they are in phase with the rest of the hem- 
isphere. The strong influence of the SH high latitudes in both modes 5 
and 6 calls for some caution in light of the sparsity of observations 
there. 

A striking aspect of the mode 6 vector is its largely uniform sign 

over the globe. In fact the sum of the 27 box values in mode 6 is larger 

than that in any other vector. This raises an important point regarding 

the relative importance of these vectors in explaining the variance in 

the global M time series. >'ur analysis approach, of course, has selected 

those eigenvectors that explain most efficiently the variance in the 

field of m D values, not in M, which is the global sum of this field, 
o 

Therefore, a vector that explains a large portion of the variance in 
the density field need not be important with regard to the behavior of 
M if that vector consists of elements that mostly balance each other, 
i.e., if the variability depicted by the vector in one part of the atmo- 
sphere is of opposite sign to equally strong variability in another part 
of the atmosphere. As we show in the Appendix, the amount of the vari- 
ance in M explained by each of our densi '» field eigenvectors is propor- 

2 th 

tional to A^p^, where A ^ is the eigenvalue of the i — vector and p^ is 

the sum of the 27 box values in the vector (which we may think of as the 

globally Integrated momentum associated with that vector). The relative 

role played by each of our six modes in terms of explaining the variance 


in M can, therefore, be readily evaluated; the result of "reordering" 
the vectors in this way is shown in Table 2. Note in this context the 
increased importance of mode 6, now just below that of mode 2. We will 
have more to say about this result in the following subsection. 

b. Variations in zonal belts 

Earlier, we mentioned our desire to compress the zonal mean wind 
data set to make our analysis more tractable. We chose to do so in the 
form of boxes because we were uncertain of the structure, both hori- 
zontal and vertical, of the EOFs that might emerge. Viewing the signif- 
icant box modes 1-6 in Fig. 13, however, it is now clear that their 
vertical structure is quite regular. This allows us to suppress the 
vertical component of variability in the data, by dealing with vertically 
integrated fields, without much loss of information. At the same time, 
we gain the freedom to use a finer horizontal resolution than in our box 
analysis while still retaining a manageable data set. Thus, we divided 
the atmosphere into as many belts as our data would allow, and com- 
puted the angular momentum of the atmosphere between 1000 and 100 mb in 
each belt for every day from 1976-80. 

Because we wished our belts to intersect equal surface areas on 
8 

the earth and because our zonal mean wind data are at 2.5° latitude 
intervals, we constrained all our belts to have the same area as the 
one between the equator and 2.5° N. ForLy-six zonal belts over the 
globe emerge in this manner. Their latitudinal boundaries are listed 

g 

The use of data that are distributed evenly over an area avoids the 
problem of "weighted double-counting" mentioned by Johnston (1981). 


27 


Table 2 


Relative importance of modes 1-6 of box analysis 
(annual signal removed) in explaining variance in M 


Mode // 

Eigenvalue (,\) 

in 48 .24-2 
10 kg m s 

M. 

X of total M 
variance explained 

2 

8.83 

1.847 

32.0 

6 

3.60 

2.852 

31.2 

4 

4.96 

1.724 

15.7 

1 

11.04 

0.957 

10.8 

5 

4.20 

0.408 

0.7 

3 

7.48 

-0.151 

0.2 

Sum 



90.6 
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in Table J. Within each belt, we computed the angular momentum of the 

atmosphere, m, via 
b 

100 mb 

lu) cos' <(> dp d4> lb) 

b 1000 mb 

where <f> runs between the southern and northern boundaries of belt b. 

As for m , we evaluated (6) by applying a trapezoidal rule (in latitude) 
n 

In such a way that ^ ra b * M each day. Also as for we deal only 
46 belts 

with 3~dav means of once-dailv m. values here. 

o 

The 5-year mean distribution of and its standard deviation are 
given in Fig. lb. As in the field of the maximum values can be 

idem it led with the zonal mean jets of each hemisphere, although we can 
now resolve the position of these maxima more precisely as lying in the 
vicinity of S and 37° N. Again, the SH maximum is larger, but the 
NH one is more variable. Most of this variance is associated with the 
annual signal, as the comparison in Fig. 17 shows. 

We performed an EOF analysis on the full nv^ time series in the same 
manner used with the box data. The first two modes (each normalized so 
that the sum of the squares of their 46 belt values is unity) and their 
associated time series are contained in Fig. 18. The first mode iden- 
tifies the spatial structure of the annual signal, and both it and mode 
2 closelv parallel their analogs for m„ in Figs. 9 and 10. 

Upon removing the annual signal, we find that the first eight modes 
pass our test of statistical significance tFlg. I $) . These modes are 
depicted In Fig, JO, their t ime series in Fig. J1 , and spectra of the 
1 1 me so r i e s in Fig. 2 2 . Mod e s 1 an d 2 c l o s e l y re semb 1 e t he i r c oun t o r - 
parts in the box analysis and little more need bo said of them. Mo Vs 3 



Table 3 


Latitudinal boundaries used to define 
46 equal area belts on the globe 


Belt 


Latitude limits (either ° 


1 

73.0 

- 90.0 

2 

65.9 

- 73.0 

3 

60.4 

- 65.9 

4 

55.7 

- 60.4 

5 

51.5 

- 55.7 

6 

47.7 

- 51.5 

7 

44.1 

- 47.7 

8 

40.7 

- 44.1 

9 

37.5 

- 40.7 

10 

34.4 

- 37.5 

11 

31.4 

- 34.4 

12 

28.6 

- 31.4 

13 

25.8 

- 28.o 

14 

23.0 

- 25.8 

15 

20.4 

- 23.0 

16 

17.7 

- 20.4 

17 

15.1 

- 17.7 

18 

12.6 

- 15.1 

19 

10.0 

- 12.6 

20 

7.5 

- 10.0 

21 

5.0 

- 7.5 

22 

2.5 

- 5.0 

23 

0.0 

- 2.5 


N or °S) 
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in the belt and box analyses are also similar, although the NH mid-lati- 
tude maximum in the latter is not so pronounced in the former, where in 
fact the largest value now occurs in the 55.7-60.4° S belt. By mode 4, 
significant differences between the belt and box analyses begin to 
appear. In particular, the maximum at the equator in box mode 4 does 
not exist in the belt mode. The fact that such differences between the 
two analyses exist should not be surprising though, in light of the addi- 
tional degrees of freedom (and number of significant modes) in the belt 
data. 

Interestingly, the major features in mode 6 of the belt analysis 
are similar to those in the box analysis mode 6. Some differences in 
detail do exist (for example in SH mid and high latitudes) , and there 
are some notable differences between the spectra of their time series. 

One important similarity belt mode 6 shares with its counterpart, though, 
is its importance in explaining non-seascnal changes in global M. In 
fact, as Table 4 shows, belt mode 6 alone explains almost half of these 
changes, more than twice the amount explained by the next largest Re- 
ordered" mode. This result is all the more remarkable when it is real- 
ized that (a) our analysis approach has not been specifically designed 
to evaluate modes of global M change and (b) mode 6 explains but 6.9% 
of the variability in m^. The implication of the latter is that non- 
seasonal changes in M on the time scales viewed here result to a consid- 
erable extent from relatively small Imbalances among regional variations 
in momentum. Much of the more important changes in momentum density, in 
the sense of large-scale cohesive changes that are depicted by modes 1-5, 
actually cancel each other in the context of a global integral. Only 
after these five modes, which in total explain 73% of the m, variance. 
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Table 4 


Relative importance of modes 1-8 of belt analysis 
(annual signal removed) in explaining variance in M 


Mode // 

Eigenvalue (X) 

, n 48 ,2 4 -2 

10 kg m s 

U 

% of total M 
variance explained 

6 

2.47 

4.186 

46.0 

2 

6.46 

1.748 

21.0 

3 

5.04 

1.704 

15.6 

8 

1.57 

-2.241 

8.4 

1 

7.29 

0-519 

2.1 

7 

1.69 

0.556 

0.6 

5 ' 

2.80 

-0.369 

0.4 

4 

4.33 

0.025 

0.0 

Sum 



94.1 
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are accounted for does the pattern emerge that explains the bulk (i.e., 

of the variability in global M (annual signal removed; Fig. 3). 

The contribution made by each belt to the non-seasonal changes in 

M can, of course, be determined directly by computing the covariance 

between the time series of and M, after removing their annual signals. 

Note that the sum of these 4b covariances will yield the variance in M 

50 2 4 - 2 

with its annual signal removed (0.94x10 kg" m s from Table 1). 

These covariances are plotted in Fig. 23, along with a plot of the 
covariance between and M when besides the annual signal, interannual 
fluctuations are also removed (by averaging the covariances determined 
separately for each of the years 1975-80). Comparing the two curves 
reveals generally small differences, indicating that most of the behav- 
ior shown is assoc in ted with non- seasonal fluctuations on time scales of 
less than a year. An exception occurs in the vicinity of 50° S where 
much, but not all, of the difference can be attributed to the 1976 anom- 
aly discussed In section 3. Finally, we note that the pattern of the 
covariance curve (annual signal removed) in Fig. 23 can be mathematically 
related to the EOF patterns in Fig. 20. Treating the covariance curve 

as a 46-component vector, we show in (All) in the Appendix that the dot 

t h 

product between it and the i — - eigenvector is given by \^u^. Referring 
to Table 4, it is apparent, therefore, that the covariance curve is 
shaped largely by modes 2, 6 and 3 (in that order). In any case. Fig. 

23 clearly illustrates that the largest contributions to non-seasonal 
variations in M on time scales of less than a year come from the sub- 
tropics of each hemisphere. The strong influence of the SH in forcing 
both interannual and higher frequency non-seasonal changes in so funda- 
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mental a parameter as M argues, of course, for increasing efforts to 
observe the atmosphere in that half of the world. 

6. Concluding Remarks 

We have diagnosed the behavior of the atmosphere's angular momentum 
during the past several years on time scales as short as a day. Our 
ability to focus on the higher frequency changes in M is a novel aspect 
of this study and was made feasible only by the development at NMC of 
routine global analyses which include the basic wind fields. Fortuit- 
ously, this development has coincided with major improvements in tech- 
niques by geodesists for measuring earth rotation. The variations we 
detected in M could therefore be checked against independent values of 
changes in the length of the day. By and large, the data sets do com- 
pare favorably on time scales of less than a year or so, thereby provid- 
ing not only a confirmation of the theory relating AM and Al.o.d. but 
also a measure of confidence in the quality of these data sets. Natur- 
ally, it would also be desirable to test the NMC -based momentum values 
against those obtained from analyses by other meteorological services, 
such as Hide et al. (1980) did for a brief period. Not only would this 
enable more definitive error bars to be placed on the calculation of M, 
but also it would provide a focus for examining differences between the 
meteorological analyses. 

We have also gone beyond the study of atmospheric momentum inte- 
grated over the whole globe to consider the manner in which fluctuations 
in momentum on regional scales are interrelated. The most obvious such 
fluctuation Involves the seasonal change in the major jets of both hem- 
ispheres, with th J asymmetry in these signals giving rise to the dominant 


34 



signature in M. To identify regional fluctuations in momentum beyonr 
those connected with the annual cycle, we adopted an empirical orthog- 
onal function technique to compute the major modes of variability in the 
momentum density field. We have not attempted here to focus in depth 
on the structure of any one particular periodicity. To do so, we might 
have considered analyzing directly the time series of the density field 
in each box or belt or, better stilJ , have pre-filtered our zonal wind 
data so as to pass only the period of interest. The latter approach has, 
in fact, been taken by Anderson and Rosen (1982) in connection with the 
^50-day periodicity. They applied a band-pass filter centered at about 
45 days to our zonal wind data and found that this periodicity appears, 
to originate in the tropics and then propagates (in phase) to NH mid- 
latitudes. (A large component of the 50-day periodicity also exists in 
the SH mid-latitude band-pass data, but unlike its NH counterpart it 
does not appear to be linked coherently with the tropical fluctuation.) 
The EOF analysis used by us here, of course, has not been designed to 
reveal such propagating phenomenon, since it utilizes only contemporary 
covariances among the momentum density field variations. 

We are inclined to look to processes within the atmosphere for an 
explanation of the various fluctuations present in the M and Al.o.d« 
time series. On the other hand, we cannot but help notice that signif- 
icant signals exist near 30 and 14 days in M and some of the EOFs of the 
density field variations, and that these periods correspond closely to 
those associated with tides. Perhaps some element of these signals in 
the atmosphere are forced by tides in the solid earth, but such specula- 
tion can be confirmed only by further detailing the atmospheric structure 
of these fluctuations. This is beyond our present scope. 
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Our last figure indicates that much of the SH accounts for a con- 
siderable portion of the non-seasonal variability in global atmospheric 
angular momentum. There has been a new appreciation among meteorologists 
for the importance of the SK in global climate, and Trenberth (1979) has 
noted the large interannual variability there. Our results for the SH, 
if generally correct, indicate that indeed a broad range of fluctuations 
in its zonal wind field exists and is important. 

Heretofore, the globally integrated value of atmospheric angular 
momentum has been more of theoretical interest to meteorologists than 
of practical consequence. It should be clear by now, however, that M 
is a fundamental measure of the interaction between the earth and the 
atmosphere. Importantly, this interaction provides an independent means 
of checking fluctuations analyzed in M. Moreover, considering the work 
by Krueger (1982) referenced earlier, it is not unlikely that M may prove 
itself worthwhile as an index of climate variability. More work needs 
to be done along these lines, of course, but we feel our efforts here 
have helped place such inquiry on a firmer foundation. 
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APPENDIX 


ORIGINAL PACC iS 
OF POOR QUALITY 


1. Variance in M explained by eigenvectors . 

Let B be the PxQ matrix of m^ (or mg) values, each of whose Q columns 
is the P-coraponent vector for a 3-day period. Thus, P ■ 27 for the box 
analysis and " * 46 for the belt analysis, and Q * 609 for either. Also, 
let u be the P-dimensional column vector all of whose elements are equal to 
1 and let M represent the Q -dimensional row vector of the values of global 
angular momentum M. Since M is the sum of all mg values, 

M - u T B (Al) 

where the superscript T indicates the transpose of a matrix (vector) . The 
variance, v , in M can be expressed as: 

V - Q -1 MM T = Q _1 u T BB T u - u T Ru (a 2) 

where R is the PxP covariance matrix formed from the mg data (i.e., R * 

q -1 bb t ) . 

Now, let E be the PxP matrix whose columns are the eigenvectors e^ of 

R, which are the EOFs derived in Section 5. Because E is an orthonormal 

T 

basis of the box vector space, then EE * I, the identity matrix, and we can 

T 

multiply the matrix R in (A2) by EE , so that 

v - u T KEE T u (A3) 

Furthermore, 

RE - EL s \ h ) 

where L is a PxP diagonal matrix whose elements are the p eigenvalues of 
R , and so 

v * U T ELE T U (A5) 
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The vector y, whose element is the sum of the components of e^ can be 
written as 


T 

y ■ E u 


0RI3I.W. . ,..•• •• 

0F POOR QJALiTY 


(A6) 


Using (A6) in (A5) we get 


v 



(A7) 


which, when written in terms of components of L and y, yields the relation- 
ship we seek: 


v 


P 

l 

i«l 


Vi 


(A8) 


2 . Relationship of eigenvectors e to the covariance between m b and M * 

The vector c, consisting of the 46 values of the covariance between m, 
- b 

and M and displayed in Figure 23 is given by 

c - Q _1 BM T - Q _1 BB T u - Ru - ELu (A9) 

T 

Multiplying both sides of (A9) by E , we get the following relation 
between two column vectors: 

E T c - Ly (A10) 

Finally, because the i^ element on the left of (A10) is simply the 
dot product of ej with g, we may write 

e^’ c ■ (All) 
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Legends for Figures 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Values of the angular momentum of the atmosphere above the 
globe, the Northern Hemisphere (NH) and the Southern Hemi- 
sphere (SH) during calendar years (a) 1976; (b) 1977; (c) 1978 
(d) 1979; (e) 1980; and (f) 1981. The scale for M is given 
along the left-hand ordinate. Inferred values of Al.o.d., 
derived from M through equation (4) , are given along the 
right-hand ordinate. Time is marked along the abscissa in 
terms of days from 1 January 1976. 

Values of the angular momentum of the atmosphere above the 
globe, as in Fig. 1 but compressed for the five-year period 
from 1 January 1976 to 31 December 1980. (Note also that the 
ordinate scale has been expanded from that in Fig. 1.) 

Three-day mean values of the angular momentum of the atmo- 
sphere above the globe during 1976-80 after removing an 
annual signal. 

Power spectrum of the angular momentum of the atmosphere above 

the globe during 1976-80 on periods of from 6 to 128 days 

2-12 

(heavy curve). The ordinate is in units of (kg ms) -day. 
The abscissa is linear with respect to frequency. The thin 
curve is an estimate of the 93% level of significance. 
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Figure 5 


Figure 6 


Figure 


Figure 


Variations ±a the length of day during 1976-80 determined 
from lunar laser ranging observations (solid line) and in- 
ferred from M through equation (4) (dashed line). The fort- 
nightly and monthly tidal terms in the lunar laser ranging 
values have been removed. The two curves are displaced 
relative to each other by an arbitrary constant amount 
(see text) . 

Variations in the length of day during 1976-79 determined 
from LLR observations (solid line) and inferred from M 
(dashed line) , as in Fig. 5 but now with the mean levels and 
annual and semiannual variations also removed each year from 
both sets of data. Time is marked in terms of days from 
1 January 1976. (Taken from Langley et al. , 1981b.) 

Mean distribution of the momentum density field m^ based on 
3-day averages of once-daily data for 1976-80. Values are 
plotted (in units of 10 Z kg in s ; at the midpoint of each 
(equal nr ss) latitude-height cross sectional box. Isolines 
have been drawn only for ease of viewing the figure; negative 
values are shaded. 

The variance of nr. in each box determined from 3-day mean 

D 

values of once-daily data for 1976-80. Units are 10^® 



Figure 9 The (a) first EOF of variability in ra D during 1976-80 and 

d 

(b) its associated time series. The percent of the total 
variance in m^ explained by this mode is given at the top. 
Values of the EOF are non-dimensional and have been multi- 
plied by 10 here; negative values are shaded. 

Figure 10 The (a) second EOF of variability in during 1976-80, 

(b) its associated time series and (c) power spectrum of 
the time series. The non-diuensional values of the EOF 
have been multiplied by 10. As in Fig. 4, the thin curve 
in the power spectrum plot is an estimate of the 95% level 
of significance. 

Figure 11 The variance of m D after removing an annual signal from 
each box’s time series of 3-day means of once-daily data 
for 1976-80. Units are 10^^ (kg m 2 s *) 2 . 

Figure 12 Percent of the variance in m^, with its annual signal re- 
moved, explained by the first 7 EOFs (solid line) and the 
same for EOFs calculated from random data (dashed line) 
for 1976-80 . 


Figure 13 EOF modes 1 through 6 [(a) through (f)] of non-seasonal 

variability in for 1976-80. The non-dimensional values 
of all the EOFs have been multiplied by 10; negative values 
are shaded. 
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Figure 14 Time series associated with each of the EOFs displayed in 
Fig. 13. 

Figure 15 Power spectrum of each of the time series in Fig. 14. The 
format and units are as in Fig. 4, with the thin curve rep- 
resenting an estimate of the 95% level of significance. 

Figure lb Mean distribution and standard deviation of the momentun 

density field derived from 3-day averages of once-daily 
data for 1976-80. The means are plotted as dots at the 
midpoint of each zonal belt. Values of ±1 standard devia- 
tion of m, in each belt are connected bv the solid line, 
b 

Figure 17 The variance in during 1976-80 derived from the full set 
of 3-day means of once-daily data (crosses) and from these 
same data but with their annual signal removed (dots). 

Figure 18 The (a) first two EOFs of variability in during 1976-80, 
(b) time series associated with mode 1, and (c) time series 
associated with mode 3. The percent of the total variance 
in m^ explained by each mode is included in its label. The 
non-dimensional values of the EOF have been multiplied by 10 

Figure 19 Percent ot the variance in m^, with its annual signal re- 
moved, explained by the t irst 10 EOFs (solid line) and the 
same for EOFs calculated from random data (dashed line) for 


1 9 76 - 80 . 


Figure 20 EOF modes 1 through 8 of non-seasonal variability in m^ for 
1976-80. The non-dimensional values of all the EOFs have 
been multiplied by 10. 

{ series associated with each of the EOFs displayed in 

4s, . 

Figure 22 Power spectrum of each of the time series in Fig. 21. The 
format and units are as in Fig. 4, with the thin curve rep- 
resenting an estimate of the 95% level of significance. 

Figure 23 The covariance between the time series of in each zonal 

belt and global M after removing their annual signals (dots) 

and after removing both their annual and interannual signals 

48 2-12 

(crosses) for 1976-80. Units are 10 (kg ms). 
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